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CHAPTER THREE 

TWO DIMENSIONAL WAVES 



1 Reflection and tranmission of sound at an inter- 
face 

Reference : L. M. Brekhovskikh and O. A. Godin: Acoustics of Layered Media I. 
Springer. §.2.2. 

The governing equation for sound in a honmogeneous fluid is given by (7.31) and 
(7.32) in Chapter One. In term of the the veloctiy potential defined by 

u = V0 (1.1) 

it is 

where c denotes the sound speed. Recall that the fluid pressure 

p = -pd<p/dt (1.3) 

also satisfies the same equation. 

1.1 Plane wave in Infinite space 

Let us first consider a plane sinusoidal wave in three dimensional space 

0(x, t) = o e l(k - x -^ = o e i(fcn ' x -^ ) (1.4) 



Here the phase function is 

0(x, t) = k ■ x — LUt 



(1.5) 
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The equation of constant phase #(x, t) = 9 Q describes a moving surface. The wave 
number vector k = kn is defined to be 

k = kn = V6 (1.6) 

hence is orthogonal to the surface of constant phase, and represens the direction of wave 
propagation. The frequency is defined to be 

—I 

Is (1.4) a solution? Let us check (1.2). 

v., ( 9 d d \ ± . w 
\ox oy oz J 



V 2 = V • V0 = ik ■ ik(f) = -k\ 



dt 2 

Hence (1.2) is satisfied if 

u) = kc. (1.8) 
Sound in an infinite space is non-dispersive. 



1.2 Two-dimensional reflection from a plane interface 

REfering to figure ??, let us cnsider two semi-infinite fluids separated by the plane 
interface along z — 0. The lower fluid is distinguished from the upper fluid by the 
subscript " 1" . The densities and sound speeds in the upper and lower fluids are p, c and 
Pi,Ci respectively. Let a plane incident wave arive from z > at the incident angle of 
9 with respect to the z axis, the sound pressure and the velocity potential are 

Pi = Pq exp[ik(x sin 9 — z cos 9)] (1.9) 

The velocity potential is 

iP 

4>i = exp [i/c(xsin$ — z cos 9\ (1-10) 

The indient wave number vector is 



= (hi, k\) = k(sm9,-cos9) 



(1.11) 
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Figure by MIT OCW. 

Figure 1: Plane wave incident towards the interface of two fluids. 

The motion is confined in the x, z plane. 

On the same (incidence) side of the interface we have the reflected wave 

p r = _Rexp[i/c(xsin6 l + zcos9)} 

where R denotes the reflection coefficient. The wavenumber vector is 

k r = (A£,A£) = k(sm9,cos9) 

The total pressure and potential are 

p = Pq {exp[ik(x sin 9 — zcos6)] + Rexp[ik(x sin9 + zcos9)}} 
iP 

(ft = {exp[ik(x sm9 — ^cos^)] + Rexp\ik(x sin 9 + zcos9)]} 

In the lower medium z < the transmitted wave has the pressure 

Pi = TP exp[iki(x sin9i — ;zcos#i)] 

where T is the transmission coefficient, and the potential 

iPo, 



p X U) 



-T exp[iki(x sin 9\ — z cos 



;i.i2) 

;i.i3) 

;i.i4) 
;i.i5) 

;i.i6) 
;i.i7) 



Along the interface 2 = 0we require the continutiy of pressure and normal velocity, 



i.e., 



p = Pi, z = 



X18) 



3.1. Reflection of sound by an interface 



4 



and 

w = w 1 , 2 = 0, (1.19) 

Applying (1.18), we get 

p Q lysine + Re ikxsme^ = T p oe ik lX sin _ < x < qo. 

Clearly we must have 

fcsinfl = jfexsin^i (1.20) 

or, 

sin 9 sin 



C Ci 

This is just Snell's law. With (1.20), we must have 



Applying (1.19), we have 



— [-Jfe cos 6e lkshie + Rk cos 9e lksine ] = — \-h cos0iTe ifelsin(?1 l 

puj p\UJ 



which implies 



Pik cos 9 

Eqs (1.22) and (1.23) can be solved to give 



pik cos 6 + pk\ cos 9 



i 



Alternatively, we have 



pccos#i + piCi cos 6* 
p\C\ cos 9 — pc cos #i 
PiCi cos 9 + pc cos #i 



Let 



Pi c 
p ' Ci 



(1.21) 



l + i? = T (1.22) 



„ pk\ cos #i 
l-i?=^ (1.23) 



T = 2pikcos6 

p/cx cos 9 1 + p! A; cos 9 

_ p 1 kcos9 - pfci cos 6>i 



^ 2piCicos9 

T = ^ (1.26) 



pici cos 6> - pc COS 01 

^ =„„„„/), „„„„„/> I 1 - 27 ) 



.'1.28) 
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where the ratio of sound speeds n is called the index of refraction. We get after using 
Snell's law that 



m cos 9 - n cos 9 1 mcos6-nJl 



sin 2 



r = 1 = y — (i 29) 

mcosfl + ncosfc mcos9+ / 1 _^e 



The transmission coefficient is 



m 2m cos 6 

T=l + R= ~^=r ( L3 °) 



We now examine the physics. 

1. If n = c/ci > 1, the wave enter from a faster to a slower medium, then 9 > Q\ and 

sin9/n < 1 always. R is real. In particular, for normal incidence 9 = 9\ = 0, 

in — n , 

R= 1.31 

m + n 

is real. If m > n > 1, i.e., pi/p > c/ci > 1, then i? is positive. If n > m > 1 then 
.R is negative. As 9 increases, 9i also increases. R decreases until at 9 = tt/2, so 
that 

#=-- = -1 (1.32) 
n 

Hence R lies on a segment of the real axis as shown in Figur l.a. or Figure l.b. 
Note that there is no reflection R = at the special incidence angle 9 = 9b, called 
the Brewster angle: 

sin9 B = -== (1.33) 
V 1 — m A 

2. If however n < 1, i.e., the wave enters from a slower medium to a faster medium, 
then 9\ > 9. for sufficiently small 9, R is real. When 9 increases to a critical value 
5, defined by 

sin<5 = n (1.34) 

sin#i = 1 so that 9\ becomes 7r/2, and 

m cos 5 



R 



m cos 5 



Below this critical angle (0 < 9 < 5), R is real. In particular, when 9 = 0, (1.31) 
applies. See Figure 2.c for m > n and in 2.d. for m < n. 
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Fig. 2.2a-d. Reflection coefficient in 
the complex plane for various rela- 
tions between the media parameters: 
G = S(*) Kc/ci < 8 i/e: 

(b) 1 < c/ci, Qilg < c/ci; 

(c) c/ci < 1, c/ci < ei/g; 

(d) qi/q < c/ci < 1 



Figure 2: Complex reflection coefficient. From Brekhovskikh and Godin §.2.2. 
When 6 > 5, the square roots above become imaginary. We must then take 

cos#i = 



sin 2 # . sin 2 6 , 
1 =- = i\l — ^ 1 



1.35) 



The reflection coefficient is now complex 

mcosfl -inyj^- - 1 



R = 



mcos^ + im/^ _ i 



;i.36) 



with \R\ = 1, implying complete reflection. As a check the transmitted wave is 
now given by 



p t = T exp 



ki (ix sin Q\ + z\j sin 2 9/n 2 — 1 



(1.37) 



so the amplitude attenuates exponentially in z as z — > — oo. Thus the wave train 
cannot penetrate much below the interface. The dependence of i? on various 
parameters is best displayed in the complex plane R = 3?i? + i^sR. It is clear 
from (1.36 ) that QR < so that R falls on the half circle in the lower half of the 
complex plane as shown in Figures 2.c and 2.d. 
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2 Equations for elastic waves 

Refs: 

Graff: Wave Motion in Elastic Solids 

Aki & Richards Quantitative Seismology, V. 1. 

Achenbach. Wave Propagation in Elastic Solids 

Let the displacement vector at a point Xj and time t be denoted by Ui(xj,t), then 
Newton's law applied to an material element of unit volume reads 

d 2 Ui drij 

P ^=d X - (2 ' 1} 
where is the stress tensor. We have neglected body force such as gravity. For a 
homogeneous and isotropic elastic solid, we have the following relation between stress 
and strain 

Tij = XckkSy + 2fieij (2.2) 
where A and /i are Lame constants and 



2 \dxj dx 

is the strain tensor. Eq. (2.2) can be inverted to give 



where 



is Young's modulus and 



'„ '//v ;^, (2.4) 



E = ^ X + ^ (2.5) 
A + /i 

A 



2(A + /i)' 

Poisson's ratio. 

Substituting (2.2) and (2.3) into (2.1) we get 

dry \ dekk 5 + d ( dui | du 3 
dxj dxj ^ dxj \dxj dxi 

de kk d 2 Ui d 2 Uj 

= \- Q \- 



(2.6) 



<9x^ dxjdxj dxidxj 

11 ■ 
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In vector form (2.1) becomes 

d 2 u 



= (A + n) V(V • u) + /iV 2 u (2.7) 



' dt 2 

Taking the divergence of (2.1) and denoting the dilatation by 



du\ du2 duz , n 

OX i OX 2 OX3 



we get the equation governing the dilatation alone 
<9 2 A 



(A + /x)V • VA + yuV A = (A + 2/i)V A (2.9) 



or, 



where 



9t 2 

«9 2 A 



9t 2 



4V 2 A (2.10) 



(2 ' n) 

Thus the dilatation propagates as a wave at the speed c^. To be explained shortly, this 
is a longitudinal waves, hence the subscript L. On the other hand, taking the curl of 
(2.7) and denoting by u the rotation vector: 

u = V x u (2.12) 

we then get the governing equation for the rotation alone 

= c 2 T V 2 Q (2.13) 



dt 2 
where 

ct = a[^ (2.14) 

V p 

Thus the rotation propagates as a wave at the slower speed ct- The subscript T indicates 
that this is a transverse wave, to be shown later. 
The ratio of two wave speeds is 



ct V A* 

Since 



Ct ' A + ">1. (2.15) 
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it follows that the speed ratio depends only on Poisson's ratio 



cl 
c t 



2-2v 
1 - 2u 



(2.17) 



There is a general theorem due to Helmholtz that any vector can be expressed as 
the sum of an irrotational vector and a solenoidal vector i.e., 



subject to the constraint that 



u = V0 + V x H 



V-H = 



(2.18) 



(2.19) 



The scalar and the vector H are called the displacement potentials. Substituting this 
into (2.7), we get 

P^[V0 + V x H] = /xV 2 [V0 + V x H] + (A + //)VV • [V0 + V x H] 
ot 



Since V • V0 = V 2 0, and V • V x H = we get 



(A + 2/i)V 2 0-p 



d 2 



dt 2 



+ V x 



/iV 2 H - p 



d 2 H 

w 



(2.20) 



Clearly the above equation is satisied if 



(A + 2/x)V 2 0-p 



and 



/iV 2 H - p 



d 2 H 

W 







(2.21) 



(2.22) 



Although the governing equations are simplified, the two potentials are usually coupled 
by boundary conditions, unless the physical domain is infinite. 
A typical seismic record is shown in Figure 2. 



3 Free waves in infinite space 

The dilatational wave equation admits a plane sinusoidal wave solution: 



0(x,t) = <f) e lk(n -*- CLt) 



(3.1) 
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Waves in Elastic Solids 




ERM 

A typical sei sinogram during an earthquake, showing 
the arrival of the P(primary), S(secondary) and R(Rayleigh) 
waves. Note that the Rayleigh waves have the largest amplitude. 



Figure by MIT OCW. 

Figure 3: A typical seismic record. P: Longitudinal wave; S: Transverse wave; R : 
Rayleigh surface wave. From Billingham & King. Wave Motion, Cambridge U Press. 
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Here the phase function is 

0(x,i) = fc(n-x-c L t) (3.2) 
which describes a moving surface. The wave number vector k = kn is defined to be 

k = kn = V0 (3.3) 

hence is orthogonal to the surface of constant phase, and represents the direction of 
wave propagation. The frequency is 

d6 

u = kc T = -— (3.4) 

A general solution is 

= 0(n • x — cit) (3.5) 
Similarly the the following sinusoidal wave is a solution to the shear wave equation; 

H = H e ife{n - X - CT * } (3.6) 

A general solution is 

H = H(n • x - c T t) (3.7) 



Note: 

We can also write (3.5) and (3.9) as 



cl 

and 



= <j>{t ) (3.8) 



H = H(t-— ) (3.9) 

c T 

where 

n n . 

s L = — , s T = — 3.10 

are called the slowness vectors of longitudinal and transverse waves respectively. 



In a dilatational wave the displacement vector is parallel to the wave number vector: 

u L = V0 = 0'n (3.11) 
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from (3.5), where 0' is the ordinary derivative of with repect to its argument. Hence 
the dilatational wave is a longitudinal (compression) wave. On the other hand in a 
rotational wave the displacement vector is perpendicular to the wave number vector, 

„ TT fdH z dH y \ (dH x dH z \ ( dH y OH, 
u T = VxH = eJ— *- )+e y — x - - — z - +e ' ' 



dy dz J \ dz dx J \ dx dy 
= e x (H' z n y - H' y n z ) + e y (H' x n z - H' z n x ) + e z (H' y n x - H' x n y ) 

= n x H' (3.12) 
from (3.7). Hence a rotational wave is a transverse (shear) wave. 

4 Elastic waves in a plane 

Refs. Graff, Achenbach, 

Aki and Richards : Quantitative Seismology, v.l 

Let us examine waves propagating in the vertical plane of x, y. All physical quantities 
are assumed to be uniform in the direction of z, hence d/dz = 0, then 

u = dl + m u = dl_dH 1 u = _dH x + dH y 
dx dy ' y dy dx dy dx 

and 
where 

^ + ^ = 1^ (4 15) 

dx 2 dy 2 c\ dt 2 ' 

PH, PH L= 1_PH L 

dx 2 + dy 2 4 dt 2 ' P X,V ' Z [ bJ 

Note that u z is also governed by (4.16). 

Note that the in-plane displacements u x ,u y depend only on and H z , and not on 
H x ,H y . Out-of-plane motion u z depends on H x ,H y but not on H z . Hence the in- 
plane displacement components independent of the out-of-plane component 
u z . The in-plane displacements (u x ,u y ) are associated with dilatation and in-plane 
shear, represented respectively by <f> and H z , which will be refered to as the P wave and 
the SV wave. The out-of-plane displacement u z is associated with H x and H y , and will 
be refered to as the SH wave. 
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From Hooke's law the stress components can be written 

f du x du y \ du x fdu x du y \ du y 



dx dy J dx \ dx dy J dy 

d 2 <p d 2 <p\ fd 2 <p d 2 H z 



dx dy J dy \ dx dy J dx 

x /<9 2 d 2 <p\ fd 2 <p d 2 H z \ . . 

T « = W^) {Txx + Tyy) = v{Txx + Tw) = x \w + w) ( } 

(du y du x \ ( n d 2 (t) d 2 H z d 2 H z \ , AnnS 



dx dy J \ dxdy dx 2 dy 2 

du z ( d 2 H x d 2 H y \ 

T *' = ^ = *\--W + djte) (421) 

*" + (4-22) 



9a; \ ddxdy dx 2 
Different physical situations arise for different boundary conditions. We shall con- 
sider first the half plane problem bounded by the plane y — 0. 



5 Reflection of elastic waves from a plane boundary 

Consider the half space y > 0, — oo < x < oo. Several types of boundary conditions 
can be prescribed on the plane boundary : (i) dynamic: the stress components only 
(the traction condition); (ii) kinematic: the displacement components only, or (iii). a 
combination of stress components and displacement components. Most difficult are 
(iv) the mixed conditions in which stresses are given over part of the boundary and 
displacements over the other. 

We consider the simplest case where the plane y = is completely free of external 
stresses, 

T yy = T xy = 0> (5.23) 

and 

r yz = (5.24) 
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It is clear that (5.23) affects the P and SV waves only, while (5.24) affects the SH 
wave only. Therefore we have two uncoupled problems each of which can be treated 
separately. 

5.1 P and SV waves 

Consider the case where only P and SV waves are present, then H x — H y — 0. Let all 
waves have wavenumber vectors inclined in the positve x direction: 

<j> = f{y)e* x - i » t , H z = h z {y)e Kx ~ wt (5.25) 
It follows from (4.15) and (4.16) that 

Z +a2f =°- W + ^^ (5 - 26) 

with 



'uj 2 /7T IT „ uj 2 



a=J-^-e = \k 2 -e, P= \h>-t 2 = \k 2 T-t 2 (5-27) 



„2 > V L S ' y \\ r- 2 



We first take the square roots to be real; the general solution to (5.26) are sinusoids, 
hence, 

On the right-hand sides the first terms are the incident waves and the second are the 
reflected waves. If the incident amplitudes A P , A s and are given, what are the properties 
of the reflected waves B P , B s 7 The wave number components can be written in the polar 
form: 

(£,«) = k L {sm6 L ,cos6 L ), {(,(3) = A>r(sin &t, cos 0t) (5.29) 

where kr) are the wavenumbers, the (6l, 0t) the directions of the P wave and SV 
wave, respectively. In terms of these we rewrite (5.28) 

_ j^^ e ik L (s\ne L x-cos9 L y-ujt) _|_ j^^ e ik L (sme L x+cose L y-u)t) ^ gg^ 

£j _ j^^^ikTisinBTX—cosOLy—Lot) _|_ q ^^ikrisin 9tx+cos Ory— 0Jt) ^ g-Q 

In order to satisfy (5.23) {r yy = r xy = 0) on y = for all x, we must insist: 

k L sin 6 L = k T sin 9 T , (f = C) (5.32) 
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cl 



c t 



implying 



sinOi cl \ + 2/i kj 



K 



(5.33) 



(5.34) 



smd T ct V n k^ 
When (5.23) are applied on y = the exponential factors cancel, and we get two 
algebraic conditions for the two unknown amplitudes of the reflected waves (Bp, Bs) : 



k 2 L (2sm 2 6 L - K 2 )(Ap + B P ) - k T sm29 T (A s - B s ) = 
k\ sin 26 L (A P - B P ) - k 2 T cos 9 T (A S + B s ) = 0. 

Using (5.34), we get 

2 sin 2 9 L - k 2 = k 2 (2 sin 2 6 T - 1) = -k 2 cos 26 t 
The two equations can be solved and the solution expressed in matrix form: 



(5.35) 
(5.36) 





\- 


Spp 


Ssp 


\B S \ 




Sps 


Sss 




(5.37) 



S = 



where 

Spp Ssp 
Sps Sss 

denotes the scattering matrrix. Thus S PS represents the reflected S'-wave due 
P wave of unit amplitude, etc. It is straightforward to verify that 

sin 26 L sin 26 T - k 2 cos 2 26 t 



Spp — 



•>PP 

Ssp 
Sps 
Sss 



sin 26 L sin 26 T + k 2 cos 2 26 t 

-2k 2 sin26» T cos 26 T 
sin 26 L sin 26 T + k, 2 cos 2 26 t 

2 sin 26 L cos 26 T 
sin 26 L sin 26 T + k 2 cos 2 26 t 

sin 26 L sin 2# T - k 2 cos 2 2^ t 
sin 2^ L sin 26 T + k, 2 cos 2 26 t 



(5.38) 
to incident 

(5.39) 

(5.40) 
(5.41) 
(5.42) 
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In view of (5.33) and 



K =c- T = Vi—Tv (5 ' 43) 
The scattering matrix is a function of Poisson's ratio and the angle of incidence. 

(i) P- wave Incidence : In this case 9l is the incidence angle. Consider the special 
case when the only incident wave is a P wave. Then A P ^ and A s = and only 
Spp and Ssp are relevant. . Note first that 9l > 9 T in general . For normal incidence, 
9 L = 0, hence 9 T = 0. We find 

Spp = -1, S PS = (5.44) 

there is no SV wave. The refelcted wave is a P wave. On the other hand if 

sin 29 L sin 29 T - k 2 cos 2 29 t = (5.45) 

then Spp = 0, hence Bp = but 13$ ^ 0; only SV wave is reflected. This is the case 
of mode conversion, whereby an incident P waves changes to a SV wave after reflection. 
The amplitude of the reflected SV wave is 

B s tan2fl r 

Ap = SPS = (5 ' 46) 



(ii) SV wave Incidence : Let Ap = but As ^ 0. In this case 6? is the incidence 
angle. Then only S$p and S$s are relevant. For normal incidence, $l = 9 T = 0, 
Sss — ~ 1) and Ssp = 0; no P wave is reflected. Mode conversion (Bp ^ 0, Bs = 0) also 
happens when (5.45) is satisfied. Since 6l > 0t, there is a critical incidence angle 9t 
beyond which the P wave cannot be reflected back into the solid and propagates only 
along the x axis. At the critical angle 

sin0 L = 1, (9 L = tt/2), or sin0 r = l/re (5.47) 

by Snell's law. Thus for v — 1/3, k — 2 and the critical incidence angle is 9t = 30°. 

Beyond the critical angle of incidence, the P waves decay exponentially away from 
the free surface. The amplitude of the SV wave is linear in y which is unphysical, 
suggesting the limitation of unbounded space assumption. 



3.5. Reflection of elastic waves from a plane boundary 



17 




Figure 4: Amplitude ratios for incident P waves for various Possion's ratios. From Graff: 
Waves in Elastic Solids. Symbols should be converted according to : A 1 — > A P , A 2 — > 
Bp, Bi — > As, B 2 — > £g. 
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Figure 5: Reflected wave amplitude ratios for incident SV waves for various Possion's 
ratios. From Graff: Waves in Elastic Solids. Symbols should be converted according to 
: A 1 — > A P , A 2 — > Bp, Bi — > As, B 2 — * B s . 
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5.2 SH wave 

Because of (5.2) 



dH x dH v n 

+ —r 1 - = 



we can introduce a stream function ip so that 



where 



Clearly the out-of-plane dispacement is 



and 



ggx + ^ = &± + ^ = v 2 ^ (5 50) 

9?/ (9a; <9x 2 <9y 2 

T ^ = % V ^ = 4dyW (5 - 51) 



The zero-stress boundary condition implies 

|=0 (, 52 , 

Thus the problem for if) is analogous to one for sound waves reflected by a solid plane. 
Again for monochromatic incident waves, the solution is easily shown to be 

ifj = [Ae-^v - Ae l(3y ) e iax ~ iut (5.53) 

where 

a 2 + [3 2 = 4 (5.54) 

We remark that when the boundary is any cylindrical surface with axis parallel to 
the z axis, the the stress-free condition reads 

r zn = 0, on B. (5.55) 

where n is the unit outward normal to B. Since in the pure SH wave problem 

zn ^ dn dn c\ dn dt 2 

Condition (5.55) implies 

^ = 0, on B. (5.56) 
on 

Thus the analogy to acoustic scattering by a hard object is true irrespective of the 
geomntry of the scatterer. 
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6 Rayleigh surface waves 

Refs. Graff, Achenbach, Fung 

In a homogeneous elastic half plane, in addition to P, SV and SH waves, another 
wave which is trapped along the surface of a half plane can also be present. Because 
most of the action is near the surface, this surface wave is of special importance to 
seismic effects on the ground surface. 

Let us start from the governing equations again 

dx 2 dy 2 c\dt 2 ' 1 ' ' 

d 2 H z d 2 H z 1 d 2 H z 

dx 2 dy 2 4 dt 2 1 ' ' 

We now seek waves propagating along the x direction 

= £ (/(y)e* x - <hrt ) , H z = ^(h{y)e l ^ tuJt ) (6.3) 

Then f(y),h(y) must satisfy 

+ K/4-c 2 )/ = °, + (^/4-« 2 )" = o, (6-4) 

To have surface waves we insist that 



a = yjt 2 ~ co 2 /c 2 L , P=yJ?- cu 2 /c 2 T (6.5) 

be real and postive. Keeping only the solutions which are bounded for y ~ oo, we get 

(f) = Ae'^e^-^, H z = Be'^e^'^ . (6.6) 

The expressions for the displacements and stresses can be found straightforwardly. 

u x = (i^Ae~ Uy -pBe~^ e 1 ^-^, (6.7) 

u y = - (aAe-^ + i^Be-^ e^ x - ut \ (6.8) 

r xx = /i j - £ 2 ~ 2a 2 ) Ae~" y - 2iP£Be-P y } e^*-^, (6.9) 

T yy = /i{(^ 2 + £ 2 ) Ae-* y + 2iPZBe-P v }e i &- ut \ (6.10) 

r xy = n\-2ia£Ae-« y + (V + f^j Be~^ e ^--*) ( 6 .n) 
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On the free surface the traction-free conditions r yy = r xy = require that 

(p 2 + e)A + 2tMB = 0, (6.12) 
-2ia£A + (ft 2 + £ 2 ) B = 0. (6.13) 

For nontrivial solutions of A, B the coefficient determinant must vanish, 

(f+e) 2 -*m 2 = o, (6.i4) 



or 

,2 



-2 



id 




(6.15) 



which is the dispersion relation between frequency u and wavenumber £. From either 
(6.12) or (6.13) we get the amplitude ratio: 



B p 2 + £2 2ia£ 
In terms of the wave velocity c = u/^, (6.15) becomes 



(6.16) 



2 -S) 2 - 4 ( i -S) I ( i -S> ' (6 ' 17) 



or, upon squaring both sides, finally 




(6.18) 



where 



c L A + 2u 2-2v 

k — 



Ct y fi V 1 — 2v 

The first solution c = u> = is at best a static problem. In fact a = (3 = £ and A = —if?, 
so that u x = u y = which is of no interest. 

We need only consider the cubic equation for c 2 . Note that the roots of the cubic 
equation depend only on Poisson's ratio, through k 2 = 2(1 — v)/{\ — 2v). There can be 
three real roots for c or cu, or one real root and two complex-conjugate roots. We rule 
out the latter because the complex roots imply either temporal damping or instability; 
neither of which is a propagating wave. When all three roots are real we must pick the 
one so that both a and (3 are real. We shall denote the speed of Rayleigh wave by cr. 
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Figure by MIT OCW. 

Figure 6: The velocity of Rayleigh surface waves cr. From Fung Foundations of Solid 
Mechanics. 



For c = 0, the factor in curley brackets is 



{.} = -16 1- 



< 



For c = ct the same factor is equal to unity and hence positive. There must be a solution 
for c such that < c < ct- Furthermore, we cannot have roots in the range c/ct > 1. 
If so, 



< 



which is not a surface wave. Thus the surface wave, if it exists, is slower than the shear 
wave. 

Numerical studies for the entire range of Poisson's ratio (0 < v < 0.5) have shown 
that there are one real and two complex conjugate roots if v > 0.263 . . . and three real 
roots if v < 0.263 .... But there is only one real root that gives the surface wave velocity 
cr. A graph of cr for all values of Poisson's ratio, due to Knopoff , is shown in Fig. 6. 
A curve-fitted expression for the Rayleigh wave velocity is 



c R /c T = (0-87 + 1 • 12i/)/(l + i/) 



(6.19) 



For rocks, A = /i and v = 4, the roots are 

(c/c T ) 2 = 4, 2 + 2/y/S, 2 - 2/Vs. 



(6.20) 
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The only acceptable root for Rayleigh wave speed cr is 

(c r /ct) 2 = (2 - 2/"s/3)3 = • 9194 



(6.21) 



or 



CR = 0.9588c T . (6.22) 
The particle displacement of a particle on the free surface is, from (6.7) and (6.8) 



u x = iA\£ 



u y = A\-a+ u 



I 3 + ^ | , n.tr-.ct, 

2£ 



ft + £ 2 \ _i(£r-u;t) 



Note that 



a = A 



g +e 2 

2£ 



> 



6 = A 



—a + 



2£ 



A 



2£ 

a-P) 2 + k 2 
2^ 



> 



hence 



and 



u x = a sm(ut — £x) , u y = b cos(ut — £x) 



a 2 b 2 



The particle trajectory is an ellipse. In complex form we have 



^ + i^L = eX p {i (cut -£x- tt/2)} 
a b 



(6.23) 



(6.24) 



(6.25) 



(6.26) 



Hence as t increases, a particle at (x, 0) traces the ellipse in the counter-clockwise direc- 
tion. See figure (6). 
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Figure by MIT OCW. 

Figure 7: Displacement of particles on the ground surface in Rayleigh surface wave From 
Fung Foundations of Solid Mechanics. 

7 Scattering of monochromatic SH waves by a cav- 
ity 

7.1 The boundary- value problem 

We consider the scattering of two-dimensional SH waves of single frequency. The time- 
dependent potential can be wirtten as 



1>(x,y,t) = 3t[<l>(x,y)e- iut ] 



where the potential <fi is governed by the Helmholtz equation 



ox z oy z ct 



(7.1) 



(7.2) 



To be specific consider the scatterer to be a finite cavity of some general geometry. On 
the stress-free boundary B the shear stress vanishes, 

/iw 2 / d(p 



T. 



-fft 



x on 



-ibjt 







hence 



dn 



0, on B 



Let the incident wave be a plane wave 

0, = Ae ticx 



(7.3) 

(7.4) 
(7.5) 
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Let the angle of incidence with respect to the positive x axis be O . In polar coordinates 
the incident wave vector is 

k = k(cos9 Q , sin# ), x = r(cos6>, sin#) (7.6) 

and the potential is 

4>i = A exp [ikr{cos 6 cos 6 + sin o sin 9)\ = Ae lkrco < e - do) (7.7) 
Let the total wave be the sum of the incident and scattered waves 

= 0/ + 05 (7.8) 

For the scattered wave the boundary condition on the cavity surface is 

<90s = _d0s Qn B / ?g x 
dn dn 1 ' 

In addition, the scattered wave must satisfy the radiation condition at infinity, i.e., it 
must propagate outward at infinty. 

Note that this boundary value problem is identical to the one for sound plane sound 
wave in air (or water) scattered by a perfectly rigid cylinder. 

We treat below the special case of a circular cavity of readius a This is one of the 
few geomentries that can be solved analytically. 



7.2 The circular cavity 

In polar coordinates the governing equation reads 

1 d ( d<f)\ | 1 d 2 | fc2 _ Q 
r dr \ dr J r 2 d8 2 

Since 0/ satisfies the preceding equation, so does <ps- 

First, It is shown in Appendix A that the plane wave can be expanded in Fourier- 
Bessel series : 

oo 

e ikrcos{e-e ) = J2e n i n J n (kr) cos n(6 - 6 ) (7.11) 

n=0 

where e n is the Jacobi symbol: 

e = 0, e n = 2, n = 1,2,3,... (7.12) 
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Each term in the series (7.11) is called a partial wave. 
By the method of separation of variables, 

<f> s (r,e) = R(r)G(6) 

we find 

r 2 R" + rR' + (fcV - n 2 )R = 0, and 0" + n 2 Q = 

where n = 0, 1, 2, . . . are eigenvalues in order that 9 is periodic in 6 with period 2n. 
For each eigenvalue n the possible solutions are 

G n = (sinn#, cos n6), 

R n ={H^(kr),H^(kr)), 

where H n 1 \kr) , H^\kr) are Hankel functions of the first and second kind, related to 
the Bessel and Weber functions by 

H^(kr) = J n {kr) + iY n (kr) } H n 2 \kr) = J n {kr) - iY n {kr) (7.13) 

The most general solution to the Helmholtz equation is 

oo 

4> s = A (A n sinnO + B n cos n6) [C n H^\kr) + D n H^(kr)} , (7.14) 

n=0 

For large radius the asymptotic form of the Hankel functions are 

Hg) ~ ./jI e *(*r-5-¥), ~ A /X e -^-f-¥) (7.15) 

V nkr n V nkr 

In conjunction with the time factor exp(— iivt), H n ^ gives an outgoing wave while H n 2) 

gives an incoming wave. To satisfy the radiation condition, we must discard all terms 

involving H n 2 K From here on we shall abbreviate H n 1 ^ simply by H n . The scattered 

wave is now 

oo 

05 = A ^ (An sin rafl + B n cos rafl) H n (kr) (7.16) 

n=0 

The expansion coefficients (A n , B n ) must be chosen to satisfy the boundary condition 

on the cavity surface 1 Once they are determined, the wave is found everywhere. In 

-'in one of the numerical solution techniques, one divides the physical region by a circle enclosing the 
cavity. Between the cavity and the circle, finite elements are used. Outside the circle, (7.16) is used. 
By constructing a suitable variational principle, finite element computation yields the nodal coefficients 
as well as the expansion coefficients. See (Chen & Mei , 1974). 
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particular in the far field, we can use the asymptotic formula to get 



4> s ~ A ^ (A n sin n6 + B n cos n6) e 



■inn/ 2 



n=0 



ikr— in/A 



irkr 



(7.17) 



Let us define the dimensionless directivity factor 

oo 

A(6) = ( A n sin nd + B n cos nO) e' inn/2 



(7.18) 



n=0 



which indicates the angular variation of the far-field amplitude, then 



<Ps ~ AA(9) 



ikr—iw/4 



Tikr 



(7.19) 



This expression exhibits clearly the asymptotic behaviour of <ps as an outgoing wave. 
By differentiation, we readily see that 



lim \fr 



kr— >oo 



Or 







(7.20) 



which is one way of stating the radiation condition for two dimensional SH waves. 
Let us complete the solution. 

Without loss of generality we can take 9 Q = 0. On the surface of the cylindrical 
cavity r = a, we impose 



d<f>i d(p s „ 

~7T + ~7T = °> r = a 
or or 



It follows that An = and 



e n i n AJ' n (ka) + B n kH' n (ka) = 0, n = 0, 1, 2, 3, . . . n 



where primes denote differentiation with respect to the argument. Hence 



J->n — Situ 1 



H' n (ka) 



The sum of incident and scattered waves is 

oo 

4> = A € n 

and 



n=0 



Ukr) - M^H n (kr) 



H' n (ka) 



cos n6 



n=0 



cos nd 



(7.21) 



(7.22) 
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ln-1 *„-5 

Figure 8: Angular distribution of scattered energy in the far field in cylindrical scattering 

The limit of long waves can be approximatedly analyzed by using the expansions for 
Bessel functions for small argument 

2 n 

J (x) ~ 1 - X - + 0{x% J n {x) ~^- r n= 1,2,3... 

2 2 n (n — IV 
Y (x)~-\ogx, Y n (x) ^——>l,n = 1,2,3... (7.23) 

7T 7TX n 

Then the scattered wave has the potential 

= \{kaf (- l -H (kr) - H^kr) cos^ + 0(kaf (7.24) 

The term H Q (kr) coresponds to a oscillating source which sends istropic waves in all 
directions. The second term is a dipole sending scattered waves mostly in forward and 
backward directions. For large kr, the angular variation is a lot more complex. The far 
field pattern for various ka is shown in fig 4. 

On the cavity surface surface, the displacement is proportional to V 2 V ; (a, 0) or 
V 2 0(a, 9). The angular variation is plotted for several ka in figure 5. 

For numerical simulations, see the website 
http: / /ocw. mit.edu/OcwWeb/Civil-and-Environmental-Engineering/l-138JFall-2004/Simulations/ 
Remark on energy flux: At any radius r the total rate of energy outflux by the 
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Figure 9: Polar distribution of 0(a, 6) on a circular cylinder. 
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scattered wave is 



I 



d9r r 



rz dt 



= /ir / d6U 



p2tt 

/ de^t 




Jo 


Or 



/jwk 4 r 
= z — 



3? [iujk 2 (f)e- iu)t ' 

do 4>* . 

or 



(7.25) 



where overline indicates time averaging over a wave period 2tt/uj. The integral can be 
evaluated by using the asymptotic expression. 

Remark: In the analogous case of plane acoustics where the sound pressure and 
radial fluid velocity are respectively, 



d(f> d(p 
P = ~ Po ~di' Ur = ~dr 



(7.26) 



the energy scattering rate is 

»2lT 



f 

Jo 



d6pU r 



u de -i 



c 



dr 



3f dB[<j> : 



c 



dr 



(7.27) 



8 The optical theorem 

For the same scatterer and the same frequency u, different angles of incidence 9j define 
different scattering problems 4>j- In particular at infinty, we have 



<j>j ~ Aj je^cosc*-*,-) + A_ e ikr-i^j (g.!) 

Let us apply Green's formula to (pi and 2 over a closed area bounded by a closed 
contour C, 

/I (*W, - fcV*,) dA = j g - ds + f c ds - is 

where n refers to the unit normal vector pointing out of S. The surface integral vanishes 
on account of the Helmholtz equation, while the line integral along the cavity surface 
vanishes by virture of the boundary condition, hence 

<90i d(p 2 



Sc is 



>2 



dn dn 



i^f )ds = (8.2) 
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By similar reasioning, we get 



ds 



d<Pl 



' dn 1 dn 
where <p\ denotes the complex conjugate of 4> 1 . 
Let us choose 0i = 2 = 4> m (8-3), and get 



ds = 



(8.3) 



(is 



c 



ds = 23 



(is ■ 



d<j>* 
dn 







Physically, across any circle the net rate of energy flux vanishes, i.e., the scattered power 
must be balanced by the incident power. 
Making use of (8.1) we get 



= 3/ rd6 



J 

Jo 



e ikrcos(e~e ) + / A (9)e ikr - i7r/4 



ixkr 



-ikcos(9 - e o ) e - ikrcos(e - d ^ - ik\ —A*J9)e- lkr+iw/4 

irkr 



= 3 / rd9 \-ikcos(9 - 9 ) + -^-(-ik)\A \ 
Jo { T^kr 

+ ^kr[co S e-e )-i] +l */±^_ ik ^ J^-A* 



+ e -ifcr[cos9-9„)-l]-W4(_ iA .) ^ _ __ Ao 

The first term in the integrand gives no contribution to the integral above because of 
periodicity. Since 3(i/) = 3(i/*), we get 







2 r Z7T 
-- / \A (9)\ 2 d9 
n Jo 

f 3 rd9 | A (-ik) + cos(^ - e o )]e t7r/4 e ikr ^' co < e - 8 ^ 



2tt 



A (9)\ 2 d9 



<^ e~ i7r/4 



d9A {9)[\ + cos(# - o )] e ikr(i-co S (e-e o )) 



For large kr the remaining integral can be found approximately by the method of sta- 
tionary phase (see Appendix B), with the result 



p2tt 

/ d&AolO) t 1 + cos (^ " 9 )}e lkr ^- co < d ~ e ^ ~ A(^ c 



2tt , 




in/ 4 



(8.5) 
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We get finally 



/ \A \ 2 d9 = -23L4o(0 o ) 



(8.6) 



Thus the total scattered energy in all directions is related to the amplitude of the 
scattered wave in the forward direction. In atomic physics, where this theorem was 
originated (by Niels Bohr), measurement of the scattering amplitude in all directions is 
not easy. This theorem suggests an econmical alternative. 

Homework For the same scatterer, consider two scattering problems 0i and (f> 2 . 
Show that 



For general elastic waves, see Mei (1978) for similar and other identities in elastody- 
namics with rigid inclusions. J . Acoust. Soc. Am. 64(5), 1514-1522. 

9 Diffraction of SH wave by a long crack - the parabolic 
approximation 

References 

Morse & Ingard, Theoretical Acoustics Series expansions. 

Born & Wolf, Principle of Optics Fourier Transform and the method of steepest de- 
scent. 

B. Noble. The Wiener-Hopf Technique. 

If the obstacle is large, there is always a shadow behind where the incident wave 
cannot penetrate deeply. The phenomenon of scattering by large obstacles is usually 
referred to as diffraction. 

Diffraction of plane incident SH waves by a long crack is identical to that of a hard 
screen in acoustics. The exact solution was due to A. Sommerfeld. We shall apply 
the boundary layer idea and give the approximate solution valid far away from the tip 
kr ^> 1 by the parabolic approximation, due to V. Fock. 

Refering to figure () let us make a crude division of the entire field into the illuminated 
zone I , dominated by the incident wave alone, the reflection zone II dominated the sum 
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Figure 10: Wave zones near a long crack 

of the incident and the reflected wave, and the shadow zone III where there is no wave. 
The boundaries of these zones are the rays touching the crack tip. According to this 
crude picture the solution is 

A Q exp(ik cos 9x + ik sin 9y), I 
A D [exp(ik cos 9x + ik sin 9y) + exp(ik cos 9x — ik sin 9y)}, II (9.1) 
0, /// 

Clearly (9.1) is inadquate because the potential cannot be discontinuous across the 
boundaries. A remedy to provide smooth transitions is needed. 

Consider the shadow boundary Ox'. Let us introduce a new cartesian coordinate 
system so that x' axis is along, while the y' axis is normal to, the shadow boundary. 
The relations between (x, y) and (#', y') are 

x = x cos 9 -\- y sin 9, y' = y cos9 — xsin9 (9.2) 

Thus the incident wave is simply 

0/ = A e lhx ' (9.3) 
Following the chain rule of differentiation, 

<90 <90 dx' <90 dy' ^ <90 . (90 
dx dx' dx dy' dx dx' dy' 
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d<f> d(p dx' d<f) dy' . d(f) ^d<p 
dy dx' dy dy' dy dx' dy' 

we can show straightforwardly that 

d 2 <p d 2 <p _ d 2 (p d 2 <p 
dx 2 dy 2 dx' 2 dy' 2 

so that the Helmholtz equation is unchanged in form in the x', y' system. 

We try to fit a boundary layer along the x' axis and expect the potential to be almost 
like a plane wave 

<l>{x: ,y') = A{x',y')e ikx ' (9.4) 

, but the amplitude is slowly modulated in both x' and y' directions. Substituting (9.4 
into the Helmholtz equation, we get 

Expecting that the characteristic scale L x of A along x' is much longer than a wavelength, 
kL x ^> 1, we have 

dA , „ , dA d 2 A 

— < kA, hence 2tk— > — 
dx dx dx z 

We get as the first approximation the Schrodinger equation 2 

, dA d 2 A , , 

m ^ + (9 ' 7) 

In this transition zone where the remaining terms are of comparable importance, hence 
the length scales must be related by 

k 1 

x' y' : 

Thus the transition zone is the interior of a parabola. 

Equation (9.7) is of the parabolic type. The boundary conditions are 



/2 -, implying ky' ~ Vkx' 



A(x,oo)=0 (9.8) 



2 In one-dimensional quantum mechanics the wave function in a potential-free field is governed by 
the Schrodinger equation 
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A(x, -oo) = A (9.9) 

The initial condition is 

am) A °; y \ >0 > (9.io) 

( A>, y' < o 

he initial-boundary value for A has no intrinsic length scales except x', y' themselves. 
Therefore the condition kL x 3> 1 means kx' 3> 1 i.e., far away from the tip. This 
problem is somwhat analogous to the problem of one-dimensional heat diffusion across 
a boundary. A convenient way of solution is the method of similarity. 
Assume the solution 

A = AJ( 1 ) (9.11) 

where 

7=^1 (9.12) 

V 7TKX 

is the similarity variable. We find upon subsitution that / satisfies the ordinary differ- 
ential equation 

/"-i7r 7 /' = (9.13) 
subject to the boundary conditions that 

/-0, 7 ^ -oo; /-l, 7-00. (9.14) 

Rewriting (9.13) as 



we get 



- = 27T7 



log /' = + constant. 

One more integration gives 



t / inn? 

f = C I exp [~Y~ ) du 



Since 



iT exp ( 



/ rcu 2 \ , e l7T / 4 
du 



2 J V2 



we get 



e -i7r/4 
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and 
/ 



A 
A, 



3 -l7r/4 ft 



cxp 



mu 



du 



in/A C p in/A 



V2 { V2 



Denning the cosine and sine Fresnel integrals by 



C( 7 ) 
we can then write 



cos 



-J7r/4 



TTV 



dv, 5(7) 



sin 



cxp 



71V 



11XU 



du \ (9.15) 



dv 



V2 



+ C(7) 



+ i 



(9.16) 



(9.17) 



In the complex plane the plot of C{^) +iS( / y) vs. 7 is the famous Cornu's spiral, shown 
in figure (9). 

The wave intensity is given by 



+ C(7) 



+ 



+ 5(7) 



(9.18) 



Since C, 5 — > -1/2 as 7 ^ —00, the wave intensity diminshes to zero gradually into 
the shadow. However, C, 5 — > 1/2 as 7 — > 00 in an oscillatory manner. The wave 
intensity oscillates while approaching to unity asymptotically. In optics this shows up 
as alternately light and dark diffraction bands. 

In more complex propagation problems, the parabolic approximation can simplify the 
numerical task in that an elliptic boundary value problem involving an infinite domain 
is reduced to an initial boundary value problem. One can use Crank-Nicholson scheme 
to march in "time", i.e., x'. 

Homework: Find by the parabolic approximation the transition solution along the 
edge of the reflection zone. 
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I Figure by MIT OCW. 

Figure 11: Cornu's spiral, a plot of the Fresnel integrals in the comiplex plane of C{pf) + 
iS(p(). Abscissa: C(j). Ordinate : iSfa). 
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Figure 12: Diffraction of a normally incident plane sound wave on a rigid half screen. 
Ordinate: f{p/) = A/A Q . Abscissa: —7 = ky' / y/irkx' . 
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10 Exact theory of Wedge Diffraction 

Refs. J. J. Stoker 1957 , Water waves, pp 120-125 
Born & Wolf, 1950 Principles of Optics 
Noble, 1963 The Wiener-Hopf Technique 

The diffraction of plane incident waves by a semi-infinite barrier is a celebrated 
problem in classical physics and was first solved by A. Sommerfeld. Several analytical 
treatments are available, including the mathematically very elegant technique of Wiener 
& Hopf is available . In the last section we gave an approximate theory by parabolic 
approximation. Here we present an exact theory for the more general case of a wedge, 
by the more elementary method of series expansion. Though the result is in an infinite 
series and not in closed form, quantitative information can be calculated quite readily 
by a computer. 

The basic ideas were described for water waves by Stoker (1957) for a wedge. Ex- 
tensive numerical results were reported by Dr. H.S. Chen (Army Corps of Engineers, 
Tech Rept: CERC 878-16, 1987). Based on Stoker's analysis more extensive numerical 
computations have been carried out by Dr. CD. Li. To facilitate the understanding 
of the physics, these results are presented in animated form in our subject website (see 
Simulations): 

http://ocw.mit.edu/OcwWeb/Civil-and-Environmental-Engineering/l-138JFall-2004/ 

We shall describe the problem for water waves; the analysis and results can be 
adapted for cracks, sound,... etc. 

Refering to Figure 10, we consider a vertical wedge of arbitrary apex angle in a 
sea of constant depth h. Let the tip of the wedge be the z axis orgin and the x axis 
coincide with one wall. The still water surface is the x, y plane. In the cylindrical polar 
coordinate system (r, 9, z) with (x = r cos 9,y = r sin 9), z, the walls are given by 9 = 
and vix with 1 < v < 2. A train of monochromatic waves is incident from infinity at the 
angle a with respect the x axis. 

In the water region defined by < 9 < un and > z > —h, the velocity potential 
$(r, 9,z,t) must satisfy the Laplace equation, 



<9 2 $ 1 <9$ 



1 d 2 $ d 2 § 



(10.1) 



dr 2 r dr 



+ 



2 002 Q z 2 



= 



r 
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x 



Figure 13: Coordinate system and subregions. 



and subject to the linearized free surface boundary conditions 

-9t 



<9$ 

~dt 



dt dz 



whicdh can be combined to give 



<9 2 $ <9$ 

dz 



dt* +9 - 



0, z = 0. 



(10.2) 
(10.3) 

(10.4) 



Along the impermeable bottom and walls, the no flux boundary conditions are 

<9$ 



dz 







9$ 

The incident wave train is given by 



on z = —h 



at 6 = and vn 



-igAo cosh k(z + h) 



r, ti e 



-ikr cos(8—a)—iujt 



oj cosh kh 

where k is the real wavenumber satisfying the dispersion relation 



(10.5) 
(10.6) 

(10.7) 



oj 2 = gk tanh kh, 



(10.8) 
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and 7r + a is the angle of incidence measured from the x axis. A is the incident wave 
amplitude. 

Because of the vertical side-walls, the three dimensional problem can be reduced to 
a two dimensional one by letting 

(10.9) 

where <p(r, 9, t) is the horizontal pattern of the velocity potential normlized for an inci- 
dent wave of unit amplitude. 

Substituting Equation (10.9) into the Laplace equation and using both the kinematic 
and dynamic boundary conditions on the free surface, the Laplace equation is then 
reduced to the Helmholtz equation 

d 2 <p 1 d<j) 1 d 2 <p , 2 „ , 

with the following boundary conditions on the rigid walls of the wedge: 

^ = at 6 = and 6 (10.11) 

o9 

The free surface displacement ( from the mean water level z=0 can be represented by 

CM,*) = ~~ = A n(r,9)e-^ (10.12) 
9 ot 

Note that n is dimensionless. 

Refering to Figure 10, the entire water region can be divided into three zones accord- 
ing to the crude picture of geometrical optics. I: the zone of incident and reflected plane 
waves, II : the zone of incident plane wave and III ; the shadow with no plane wave. 
This crude picture is discontinuous at the border lines separating the zones. Siknce the 
physical solution must be smooth every where we must find the transitions. Let us use 
the ideas of boundary layers. 

Let the total potential be expressed in a compact form by 

= <f) (r,6) +<j) a (r,8), for all < 9 < vix (10.13) 

where (p consists of only the plane waves, 

4>i + 4> r it — a > 9 > 0, in I; 
(f) (r,9) = { fa n + a>9>n-a, in II; (10.14) 

9 > 9 > 7i + a, in III. 
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Here fa denotes the incident wave 



fa = e 



-ikr cos(6—a) 



where a denotes the angle of incidence, and fa the reflected wave 



4> r = e 



—ikr cos(8+a) 



(10.15) 



(10.16) 



The correction is the diffracted wave <p s which must satisfy the radiation condition and 
behaves as an outgoing wave at infinity, i.e., 



or 



lim \fr(^- - ikfa) = 
r^o or 

A(6)e ikr 

at r — > oo 



kr 



(10.17) 



(10.18) 



10.1 Solution by Fourier series (or, finite Fourier Transform) 

Let us solve the scattered wave formally by Fourier series 



1 2 °° 

far, 6) = —fo( r ) + — J2$ n 

U71 Z/7T ^ 



nd 



r cos 



10.19) 



n=l 



then, the Fourier coefficients are: 

4> n (kr) 



nd 

fakr, 9) cos — d9 



v 



10.20) 



From (10.10), each Fourier coefficient satisfies 



,d 2 <p n d<p n 

h r h 

dr 2 dr 



w-G)' 



The general solution finite at the origin is 



(p n (kr) = a n J n/y (kr) 

where the coefficient's a n ,n = 0, 1, 2, 3, ... are to be determined. 
The Fourier coefficient of (10.13) reads 

f un n8 f un n6 
4>n(kr) = a n J n / u (kr) = / S cos — d6 + / o cos — d6 

Jo v Jo v 



(10.21) 



10.22) 



(10.23) 
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or 



4> s ,n = a n J n/u (kr) - (p 0tn (10.24) 

Applying the operator lim^oo y/r(d/dr — ik) to both sides of (10.23), and using the 
Sommerfeld radiation condition (10.17), we have 



lim \fr ( 

r-+oo \ 



f d_ 

dr 



— ik 



a nJn/u(kr) 



, cos — dU 







10.25) 



We now perform some asymptotic analysis for large kr to evaluate a Ti 
The first part can be treated explicitly for large kr, since 



J n /u{kr) 



7i kr 



riTT 7T\ 

cos kr 

2v 4/ 



It follows that 



lim \fr 



d_ 
dr 



ik ) J n/u (kr) 




2/c e _ i ( fcr _- +f ) 

7T 



10.26) 



10.27) 



For the second part, we substitute 4> from (10.15) and (10.16) to rewrite the integral as 

1 2 



,cos — d6 



v 



-J 

Jo 



e -ikr COS (e- a)cos ™ dd + 

V 



e- lfercos(e+Q) cos— d9 
v 



f-K+a n 

+ / e - ikrcoB ^cos— d6 

v 



10.28) 



Each of the integrals above can be evaluated for large kr by the method of stationary 
phase (again). Details are given in the next subsection, only the results are cited below. 
The first integral is appoximately 



h {B) = cos (™) e~ ik ^ 



2tt' 
kr 



from which 



lim \f i ( — ik 



kr 



nO 



10.29) 



e -ikrcos(e-a) cos _ d Q 
V 



lim \fr 



d_ 
dr 



— ik I < cos 





"2tt" 


n 









2v / 2^cos(^)e- lfcr - 



(10.30) 
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where we have used i = e m l 2 . By similar analysis the second integral is found to be 



h{9) ~ 7. cos 



It follows that 



lim ^/r 



d_ 
dr 



— ik 



lim \fr 



d_ 
dr 



ik 



2ir 
kr 



v 

lcos( n(7r ; a V - j 



2tt 
kr 



= 



(10.31) 



10.32) 



Finally the third integral is approximately 



hence 



1 + a )_ykr-f 



J 3 (0)w-cos( 

_ IS 



2tt 
kr 



lim 



d_ 
dr 



ik 



V 



lim \fr { — ik 



\dr 



1 n(ir + a) ifcr _ 
-cos( — -)e* fer 

2 v 1/ ; 



in 


"2tt" 


n 


4 


_ kr 









(10.33) 



(10.34) 



In summary, only the first integral associated with the incident wave furnishes a 
nonvanishing contribution to the expansion coefficients, i.e., 

lim ( - ik) r <Po cos —d9 ~ 2v / 27fc cos !^ e ~<^+f ) (10.35) 

r-+oo \Or J J Q V V 

With this result we get by substituting (10.27) and (10.35) into (10.25), the coefficients 
a n are found 

a n = 2vr cos — e ^ (10.36) 
By the inverse transform, (10.19), we get the exact solution, 

2 



<j>{r,0) 



Jo(kr) + 2 e 2 " J n /u\kr) cos — cos 

n=l ^ 



V 



(10.37) 



Numerical computations by MATLAB is straightforward. Aside from its own practical 
interest, this solution is useful for checking strictly numerical methods for diffraction 
problem for other geometries. 
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10.2 Asymptoic approximation of Integrals 

For the first integral ii, we take the phase to be fi(0) = kcos(9 — a). The points of 
stationary phase must be found from 

f[(0) = -ksm(9-a) = 0, (10.1) 

hence 9 = a, a ± it. Only the first at 9± — a lies in the range of integration (0, n — a) 
and is the stationary point. Since 

ft fa) = —k cos(9 1 - a) = -k < (10.2) 

the integral is approximately 



h (e) w cos(^)e-^ cos ^- a ) + ^ 



2tt 



cos( — )e 



-ikr+- 



2ir 
kr 



(10.3) 



For the second integral I 2 , we take the phase to be /2(d) = kcos(9 + a). The 
stationary phase point must be the root of 

f' 2 {Q) = -ksm(6 + a) = (10.4) 

or 6 = —a, —a ± 1. The stationary point is at 62 = vr — a which is the the upper limit 
of integration. Since 

f»(9 2 ) = -k cos(e 2 + a) = k > (10.5) 

I 2 is approximately 



7 2 (^) Ri Icos(^)e-^ cos ^ +a )-^ 



2tt 
kr 



1 nin — a) 
- cos( 

2 v 1/ 



jfcr- 



2tt 
/cr 



(10.6) 



Lastly for the third integral ^3, the phase is fz{0) = kcos(9 — a). The point of 
stationary phase is found from 

f' 3 (9) = -kwi(e-a)=Q (10.7) 

or 9 = ii, ±7T + a. Only the point 9 3 = n + a is acceptable and coincides with the upper 
limit of integration. Since 

ft(9 3 ) = —k cos(0 3 — a) — k > 0, (10.8) 

I 3 is approximately 

l ^<<* + %^\^ (1Q9) 



/ 3(^)~^C0S( 
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10.3 Two limiting cases 

(1) A thin barrier. Let the wedge angle be by setting v — 2. Equation (10.37) then 
becomes 

oo „ 

0(r, 9) = Jo(kr) + 2 £ e'^ J n/2 (kr) cos ™ cos ™ (10.10) 

n=l 

(see Stoker (1957)). 

(2) An infinite wall extending from x = — oo to oo. Water occupying only the half 
plane of y > and the wedge angle is 180 degrees. The diffracted wave is absent from 
the solution, and the total wave is only the sum of the incident and reflected waves : 

0( r , 6) = e - ifcrcos ( e - Q ) + e -i*rc«.(0+a) (10.11) 

By employing the partial- wave expansion theorm, (Abramowitz and Stegun 1964), the 
preceding equaion becomes 



4>(r,e) = 2 



;io.i2) 



Jo(kr) + 2 ^^(— i) n Jn(kr) cosnacosnO 

n=l 

which agree with (10.37) for v — 1. 

For the animated version of the wave patterns for different wedge angles snd angle 
of incidence, please visit the website. 
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A Partial wave expansion 



A useful result in wave theory is the expansion of the plane wave in a Fourier series 
of the polar angle 9. In polar coordinates the spatial factor of a plane wave of unit 
amplitude is 



Consider the following product of exponential functions 



e zt/2 e -z/2t 



n\ V 2 



n=0 



n\ V It 



£ 1 "'' 

n=0 

*/2)» (,/2f+ 2 + + + ( 1)r (z /2)"+* + 



n! l!(n + l)! 2!(n + 2)! 
The coefficient of t n is nothing but J n (z), hence 

1 



r\(n + r)\ 



cxp 



Now we set 



t = ie z = kr. 



The plane wave then becomes 

oo 

e ikx = e m{e+7T/2 J n {z). 

N=-oo 

Using the fact that J_ n = (— l) n J n , we finally get 

oo 

e ikx = eJ fcrco S e = Y^e n i n J n {kr) cosn6, 



n=0 



(A.l) 



where e n is the Jacobi symbol. The above result may be viewed as the Fourier expan- 
sion of the plane wave with Bessel functions being the expansion coefficients. In wave 
propagation theories, each term in the series represents a distinct angular variation and 
is called a partial wave. 

Using the orthogonality of cosnO, we may evaluate the Fourier coefficient 



Jn(kr) = — 

€ n l""K 



re"™ 7 ' JO 



Akr cos 8 



cos n6d6, 



(A.2) 



which is one of a host of integral representations of Bessel functions. 
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B Approximate evaluation of an integral 

Consider the integral 

/•27T 

/ dd[l + cos(# - o )] e ^r(i-cos(e-e o )) 
Jo 

For large kr the stationary phase points are found from 

^[l-cos(0-0 o )]=sin(0-0 o ) = O 

or 9 = 9 Q , 9 D + tc within the range [0, 2tt}. Near the first stationary point the integrand 
is dominated by 

2A{9 )e^ d - d ^' 2 . 

When the limits are approximated by (— oo, oo), the inegral can be evaluated to give 

A(9 ) e ikr62 ' 2 d9 = x [^e^ 4 A(9 ) 
J -oo V kr 

Near the second stationary point the integral vanishes since 1 + cos(# — 9 ) — 1 — 1 = 0. 
Hence the result (8.5) follows. 



